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Abstract 

Flash Analog to Digital Converter (FADC) is frequently used in nuclear and 
particle physics experiments, often as the major component in big multi-channel 
systems. The large data volume makes the optimization of operating parameters 
necessary. This article reports a study of a method to extend the dynamic range of 
an 8-bit FADC from the nominal 2*^ value. By comparing the integrated pulse area 
with that of a reference profile, good energy reconstruction and event identification 
can be achieved on saturated events from CsI(Tl) crystal scintillators. The effective 
dynamic range can be extended by at least 4 more bits. The algorithm is generic 
and is expected to be applicable to other detector systems with FADC readout. 

PACS Codes: 07.05.Kf, 07.05.Pj, 07.50. Qx 
Keywords: Data analysis. Image Processing, Electronics 



( Submitted to Nucl. Instrum. Methods A ) 



* Corresponding author: Email: htwong@phys.sinica.edu.tw; Tel:+886-2-2789-6789; FAX:+886-2-2788- 
9828. 



1 Introduction 



With the advance of electronics, data acquisition (DAQ) and computing capabihties, it 
is more and more common to use Flash Analog to Digital Convertors (FADCs) as the 
general-purpose readout device for different detectors. Typically, the entire pulse shape 
from the individual detector module, as well as the relative timing between pulses from 
different detector channels, are recorded. Such readout scheme replaces Analog Digital 
Convertors and Time-to-Digital Convertors, which give information only on the amplitude 
or integrated area and the timing of the pulses, respectively. 

With FADCs as readout, the information related to a single pulse is represented by 
many measurements. Frequently, multi-channel (order of hundreds and thousands is com- 
mon) FADC systems are deployed. As a result, a large volume of data have to be ex- 
tracted, stored and processed. The n-bit resolution per time-bin and the sampling rate of 
the FADC are important parameters which will have tremendous effect on the character- 
ization and performance of the experiment such as detector resolution, DAQ dead time, 
data size, data processing time, cost and so on. 

A motivated question to address is therefore 'What is the dynamic range of an n-bit 
FADC, beyond the nominal 2^^ value?" That is, can the crucial information like the total 
area as well as the rise and fall times be reconstructed, even though part of the pulse — 
the largest and usually the most important part — is not recorded? These are the focuses 
of investigation of this article. 

In Sectional we describe the research program on which this project is based, and the 
hardware used to generate the data which allow such an investigation. Two approaches 
are surveyed and compared. One of the algorithms, the Charge Summation Correction 
Method, showed better potentials and were pursued in depth. The goals are to achieve 
both energy reconstruction and particle identification even for saturated events. The 
details and results are discussed in Sectional 

Though the algorithm was developed and tested with scintillating Csl(Tl) crystal 
detectors, the approach is generic and can be applicable to other detector schemes using 
FADC readout as well. 

2 Scintillating CsI(Tl) Crystal Detector 

The present investigation is part of the research program of the TEXONO Collabora- 
tion pi. A CsI(Tl) scintillating crystal [2j array with 2 kg modular mass and a total 



mass of 186 kg is taking data at the Kuo-Sheng (KS) Nuclear Power Plant to study low 
energy neutrino physics. Each module is read out by two photo-multipliers (PMTs) and 
the signals are recorded by a 20 MHz FADCs with 8-bit resolution The FADC value 
at the t[^ time-bin is denoted by yi. Throughout this work, we use the normalized value 
Yj which is yi divided by the saturation level Lsat = 2" — 1 = 255 : Yj = yi/Lsat • 

The physics requirements of the experiment are that events with energy from 10 keV to 
10 MeV can be reconstructed, and that pulse shape discrimination (PSD) can be applied 
to these events to achieve event (a/7) identification. That is, a dynamic range of at least 
1000 is necessary with an 8-bit resolution readout per time-bin. 

Special data were taken to devise techniques and to develop software algorithms ad- 
dressing this issue. A 5 x 5 x 5 cm^ CsI(Tl) crystal was used and signals were read out 
by a single PMT and a FADC channel at identical settings to the KS experiment. The 
trigger was set at the to = 100**^ FADC time-bin, such that the pulse information of 
5t_=5 fis (100 time-bins) and 5t+=25.6 /is (512 time-bins) before and after the trigger 
were recorded, respectively. Both a and 7 events were taken by activating the CsI(Tl) 
crystal with ^^^Am-Be and ^^"^Cs sources, respectively. 

The amplitude of the pulse was adjusted by selecting the appropriate PMT high 
voltage. Typical unsaturated and saturated pulses are displayed in Figures and [T)d, 
respectively. A large sample of unsaturated events are added together to produce a smooth 
reference pulse shape Y(t), where the peak amplitude is normalized to 1. The profiles of 
Y(t) due to a- and 7-events are depicted in Figure |2l and can be fitted to the function P] 

Y(t) = N * [ 1 - exp(— ^) ]*[- exp(--) + - exp(--) ] . (1) 

To Ti Ti T2 T2 

The rise time(ro) and fall times(ri, T2) as well as the ratio between the slow and fast decay 
components(r) are listed in Tabled The normalization constant N is chosen to make the 
peak amplitude equal to unity. The values of tq are dominated by the electronics shaping 
rise time of 250 ns [3], while (ri, T2) are due to the characteristic fluorescence time profile 
of CsI(Tl). 

The integrated charge (Q) of an event is obtained by numerically integrating the 
pedestal-subtracted FADC pulses: 

Q = E [Y,-Yo] = E , (2) 

ti>to ti>to -"^sat 

where yo is the average pedestal level determined from the pre-trigger part (ti < to) of 
the pulse, and Yq = yo/Lsat- The typical full-width-half-maximum (FWHM) resolution 
for unsaturated events at 662 keV is about 8%. 



3 Charge Summation Method 



3.1 Initial Survey 

Saturated events are those with pulse amphtude larger than those corresponding to the 
maximum FADC value of Lgat = 255. The focus of this section is to devise an algorithm 
with which energy reconstruction and event identification by PSD can be achieved for 
these class of events. 

Two approaches were surveyed as initial studies: 

1. Statistical Approach In principle, one can derive the values of N in Equation d 
using maximum likelihood analysis or minimum fitting methods, based only on 
the (ti,yi) data where yi < Lsat- The values of (ro,ri,r2,r) are adopted from those 
in Table ^ In practice, the fits are computer-time consuming and unstable. The 
efficiencies are bad, especially when the saturation is severe. The main reason is that 
the values of tq and ri are derived from small portions of the pulse shape. When 
the data from these portions are missing due to saturation, there will be insufficient 
constraints such that it is difficult for the analysis procedures to converge to the 
optimal values for N. 

2. Pulse Analysis Approach The objective of this approach is to identify measure- 
able values which are related to the physics variables of interest, namely, the total 
energy and the decay time. The observables are the "saturated charge" and "mean 
time". This method uses all the available information from the saturated pulses. 
Such an approach can be compared to — and in fact inspired by — the "Double 
Charge Method" |3] commonly-used to achieve a/7 identification in liquid [Sj and 
crystal j2] scintillators. 

Results of the initial surveys indicated that the pulse analysis approach would show 
better promises. The "Charge Summation Method" based on this approach was therefore 
pursued in details. We discuss the algorithm, its consistency checks and performance in 
energy reconstruction and event identification in the following sub-sections. 

3.2 Reference Profile and Correction Function 

From a standard reference pulse shape, a relationship is obtained between the "saturated 
charge (Qsat)" defined by the integration of the saturated events and the "total charge 
(Qtot)" from that of the raw fully-recorded pulses. 



The saturation level for an n=8-bit FADC is Lsat = 2° — 1 = 255. To investigate the 
saturation effects on a pulse shape with amphtude A > Lgat, a profile Z(A) is constructed 
from the reference shape Y(t) of Figure |2l scaled-up by a factor of of A/Lgat- Therefore, 
the value at the t^'^ time-bin is given by 

Zi(A) = [ A ] . (3) 

^sat 

The pedestal level Zq is the average of Z; for tj < to- The dependence of the total and 
saturated charge as functions of A for 7-events can be defined as: 

Qtot (A) = E [Zi-Zo] (4) 

ti 

and 

6 = / for z!>l ' 

respectively. That is, given the profile Z(A) with peak amplitude A/Lsat, Qtot and Qsat 
are the integrated areas of the complete pulse and the pulse truncated at a maximum 
amplitude of A = Lgat, respectively. By varying the values of A/Lgat, the dependence of 
Qtot and Qsat on A/Lgat for 7- and a-events can be derived. They are depicted in FigureEK 
and Eb, respectively. As expected, the variations shows a linear relationship between Qtot 
and A/Lgat, while Qgat would level off at large A/Lgat due to saturation. 

A relationship between Qtot and Qgat for the saturated events (that is, when Qtot > 
Qsat) can then be obtained. It is described by a polynomial function 

3 

Qtot = E '^k Qsat^ (6) 
k=0 

for both the 7- and a-samples as displayed in Figured The coefficients ai are obtained 
from minimum x^-fits. 

This function is then used to compensate the missing area in saturated events. Inte- 
gration of the saturated pulses provide the values of Qsat, from which the correct values 
of the total charge Qtot can be derived in an event-by-event basis. 

3.3 Consistency Tests 

The correction algorithm discussed in the previous section does not depend on the exact 
saturation level Lgat. One can use this feature to perform consistency checks on data 
sample where the total charge is known. An unsaturated event is made saturated by 
software. The algorithm is applied to the simulated saturated event, and the corrected 
total charge is compared with the raw measurements. 



A level of Ltest=100 is chosen. A new set of is obtained by going through Eqs. El 
to El again, with Lgat replaced by Ltest- A data set of 7-events was taken with a ^^'^Cs 
source. The PMT high- voltage was adjusted such that the amplitude for the 662 keV line 
is below the actual saturation level of 255 but significantly above Ltest- The charge Q is 
measured by integrating the normalized FADC counts Yi, and the spectrum is displayed 
in Figure E^. The same data is then made saturated by re-defining Y; Y[ such that 

Y/ _ I Yi for Yi < Ltest /Lsat 

[ Ltest /Lgat for Yi > Ltest /Lsat 

The saturated charge is derived by integrating the normalized Y( and the spectra is shown 
in FigureEb, while that for the corrected charge Q' after the correction procedures is shown 
in Figure Et- There is excellent agreement on the ^^^Cs peaks in terms of their integrated 
charge, resolution and area between Figures and|St. The event-by-event deviations are 
parametrized by 

^ - (Q + Q')/2 ■ 
The distribution of A is depicted in Figure IHl The results show that the mean bias is less 
than 0.5% which is small compared to the intrinsic FWHM resolution of 8.7%. Therefore, 
the correction procedures are valid and justified. 



3.4 Results 

3.4.1 Energy Reconstruction 

The algorithms are applied to actual saturated data. The data is generated by increasing 
the PMT high voltages such that the events corresponding to the 662 keV 7-line from 
^^''Cs have amplitude larger than Lsat=255, while those from the 32 keV line remain 
unsaturated. 

An example of the energy spectra of the saturated events, as well as the corrected 
spectra are illustrated in Figure jT^i and Figure |7|d, respectively. The ratio of the peak 
positions at 32 keV and 662 keV (denoted by Q32 and Q662, respectively) is measured and 
compared to the nominal values: 

^ 032/0662 _ ^ ^ 

32/662 ^ ^ 

This can be used as a figure of merit to quantify the performance. The values of R for 
saturated and corrected events as functions of the expected peak amplitudes are displayed 
in Figure |HK- It can be seen that the charge summation correction methods extend the 



effective dynamic range of the FADC by at least a factor of 14, or equivalently, 4 more 
bits. 

The dependence of performance of the algorithm on the FADC recording period 5t+ is 
studied. The data are truncated by software to a shorter interval and the same correction 
procedures are repeated. The results are displayed in Figure HJo for the variation of R with 
peak amplitudes at 5t+ = 25.6, 12.8, 6.4 fis (that is, 512, 256 and 128 time-bins, respec- 
tively). It can be seen that the charge reconstruction is no longer valid for 5t+ = 6.4 /is at 
an amplitude of >8 times of Lsat- Alternatively, a plot of R versus 5tunsat/'^t_,_ is displayed 
in Figureini where (5tunsat is the unsaturated interval, where the FADC counts are less than 
Lsat- It can be seen the values of R deviate from 1 only for 5tunsat/<5t+ < 0.02. Accord- 
ingly, the algorithm is valid so long as the the unsaturated interval of the event is more 
than 2% of the FADC recording period. In comparison, the fitting method as discussed in 
Section ITTl would start to fail when the saturation interval is more than twice the value 
of Ti, or equivalently, the figure-of-merit R deviates from 1 when 5tunsat/^t_,_ < 0.74. 

The dependence of energy resolution of the 662 keV line on amplitudes is shown in 
Figure IHt- There is no intrinsic bias in the correction, though the energy resolution gets 
worse from around 9% to 16% FWHM when the amplitude increases by a factor of 14 
at 5t+=25.6 fis. The degradation in energy resolution is larger for the shorter FADC 
recording time. 



3.4.2 Pulse Shape Discrimination 

As depicted in Figure|21 the scintillation light in CsI(Tl) has different decay time constants 
for 7- and a-activated events. The a-events decay faster in inorganic crystal scintillators 
— opposite to the responses in organic liquid scintillators OEj. Accordingly, the pulse 
shape information recorded by the FADC can provide event identification capabilities. 

Identification of 7/a events is achieved by studying and comparing the different re- 
sponses from two observables. One of them, the saturated charge (Qsat), has been adopted 
in energy re-construction as discussed in Section |21 The other is the "mean time" (< t >) 
which reflects the decay profiles, defined as: 

E [ti-to]*[Yi-Yo] 

E |Y,-Y„) ■ (^°) 

ti>to 

The average variations of < t > with Qsat for both 7- and a-events can be derived from 
the reference pulse shape of Figure El and are displayed in Figure ITUk . The scattered plot 
of < t > versus Qsat for individual events is shown in Figure ITHb . It can be seen that 



a/7 events are well separated for Qsat < 120. From Figure El it can be concluded that 
0/7 differentiation can be achieved for a-events with amplitude up to 20 times of the 
saturation level. 

4 Summary 

An algorithm to process saturated events in Flash Analog-to-Digital Convertor was devel- 
oped and studied. Charge reconstruction and event identification via pulse shape analysis 
are achieved, effectively extending an 8-bit FADC by at least 4 more bits. The algorithm 
is valid so long as the unsaturated interval is more than 2% of the FADC recording period. 

Although this work originates from the goals of handling data from CsI(Tl) crystal 
scintillators, the approach should be generic and can be applicable to other detector 
systems where the entire pulse shape are recorded while part of which are saturated. These 
studies will facilitate the specification and design of large FADC-based data acquisition 
systems. 

This work was supported by contracts 90-2112- M-001-037 and 91-2112-M-001-036 from 
the National Science Council, Taiwan, and 19975050 from the National Science Founda- 
tion, China. 
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Event Type 


Rise Time 


Decay Time Constant 


Ratio (r) 




[^0 (ns)] 


Fast Comp. 


Slow Comp. 

[7-2 (ms)] 




a 


203±3 


0.54±0.10 


2.02±0.02 


0.29±0.02 


7 


261±2 


0.87±0.10 


5.20±0.04 


0.61±0.01 



Table 1: Fitted rise and decay time constants as well as the ratio between slow and fast 
decay components for a and 7 events measured by CsI(Tl) crystal. 
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Figure 1: Typical (a) unsaturated and (b) saturated pulses recorded by the FADC. 
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Figure 2: The reference pulse shape for CsI(Tl) from unsaturated events due to 7-rays 
at 662 keV and a-particles at 5.4 MeV. 
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Figure 3: The dependence of (a) the total charge Qtot and (b) the saturated charge 
on A/Lsat- 




Figure 4: The total charge Qtot as a function of the saturated charge Qsat, for both 7 
and a events. The data are fitted to a polynomial function for parametrization. 




Figure 5: Energy Spectra for ^^'''Cs: (a) Raw spectra of unsaturated events, (b) Spec- 
tra of events made saturated by software, and (c) Spectra of saturated events after the 
corrections are applied. 




Figure 6: The deviation of the charge measurements between the raw unsaturated events 
and the saturated events with corrections apphed. 




Figure 7: The energy spectra of saturated events (a) before and (b) after the correction 
algorithms are apphed. 
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Figure 8: (a) The ratios of peak positions between the 662 keV and 32 keV 7-hnes before 
and after saturated pulse correction; (b) the ratios of peak positions at different time- 
bin records 5t+, and (c) the energy resolution at 662 keV for the corrected spectra, as 
functions of the expected peak amplitudes. 




Figure 9: Variation of the normalized peak ratio R with (5tunsat/<^t+, where St+—6A /is 
and 5tunsat is the mean unsaturated interval for the events. 




Figure 10: The relationship between < t > and Qsat for 7- and a-events: (a) the average 
value and (b) event-by-event scattered plot. 



